Observation of Stochastic Coherence in Coupled Map Lattices 
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Chaotic evolution of structures in Coupled map lattice driven by identical noise on each site 
is studied (a structure is a group of neighbouring lattice-sites for whom values of dynamical vari- 
able follow certain predefined pattern). Number of structures is seen to follow a power-law decay 
with length of the structure for a given noise-strength. An interesting phenomenon, which we call 
stochastic coherence, is reported in which a rise of bell-shaped type in abundance and lifetimes of 
these structures is observed for a range of noise-strength values. 



oo , 
0^ ■ 

On ■ 

1—5 

in 

(N 

> ■ 
\o ■ 

o ; 
o 

(N 
O 

\o . 

"a : 
o 



X 



Coupled map lattice (CML) model has been observed 
to exhibit diverse spatio-temporal patterns and struc- 
tures Q, and serves as a simple model for experimen- 
tal systems such as Rayleigh-Benard convection, Taylor- 
Couette flow, B-Z reaction etc. Q. A major interest is in 
understanding formation of structures, localised in both 
space and time, in turbulent fluid Role of fluctuations 
in onset, selection and evolution of such patterns and 
structures has been studied in some detail j|J^. In this 
communication we report a novel phenomenon observed 
in the dynamics of structures in a chaotically evolving 
one dimensional CML driven by identical noise. By a 
structure we mean a group of neighbouring sites whose 
variable-values follow certain prespecified spatial pat- 
tern. Distribution of these structures during evolution 
of the lattice shows that their numbers exhibit power- 
law decay with length of the structure for a given noise- 
strength, with an exponent which is a function of noise- 
strength. It is observed that average length of these 
structures shows a bell-shaped curve with a character- 
istic peak, as a function of noise-strength. Similar be- 
haviour is observed for average lifetime of these struc- 
tures during their evolution, within same range of noise 
values. We call this phenomenon stochastic coherence. 

We consider a one dimensional CML of the form 
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where Xt(i), i = 1, 2, • • ■ , L, is value of the variable lo- 
cated at site i at time t, rjt is the additive noise, s is 
the (nearest neighbour) coupling strength, and L is the 
size of the lattice. Logistic function F{x) = fix{l — x) is 
used as local dynamics governing nonlinear site-evolution 
with nonlinearity parameter /i. We have used both open- 
boundary conditions xt{0) = xt{l), Xt{L -I- 1) = xt{L) 
, and periodic-boundary conditions Xt{L + i) — Xt{i), 
for our system. For noise r]t we have used a uniformly 
distributed random number bounded between —W and 
+W, where W is defined as noise-strength parameter. 

We define a .structure as a region of space in which 
the dynamical variables at sites within this region follow 
a predefined spatial pattern M. To study coherence in 



the system we choose a spatial pattern where the differ- 
ence in the values of the variables of neighbouring sites 
within the structure is less than a predefined small pos- 
itive number say S, i.e., \xt{i) — Xt{i ± 1)| < 5. We call 
5 the structure parameter. We look for such patterns, or 
coherent structures, to appear in course of evolution of 
the model given by Eq. |l|. 

We now present our main observations. Values of /i, 
£ and L are chosen so that the resultant dynamics of 
the system is chaotic. Coherent structures with length 
< 3 (sites) and lifetime < 2 (timesteps) are disregarded 
in our observation. Values for W are chosen within the 
range [0,1]. Figure 1 shows plot (on log-log scale) of 
distribution of number n{l) vs. length I of structures for 
different values of W, with = 4, e = 0.6, 5 = 0.0001 
and L — 10000, and open-boundary conditions are used. 
Power-law nature of decay of n{l) is clearly evident, 
which has a form 



n{l) oc I 



(2) 



where a is power-law exponent. This indicates that the 
system does not have any intrinsic length scale. It may 
be noted that in absence of noise {W — 0) the decay is 
manifestly exponential Exponent a (i.e., slope of the 
log-log plot in Fig. 1) is seen to depend on noise-strength 
W . This fact is corroborated in Fig. 2 which shows plot 
of a vs. W for the same parameter values as in Fig. 1. 
The exponent is exhibiting a clear minimum for values 
of W around 0.6. We define average length I of a struc- 
ture as I = '^Inil)/ ^n{l). In Fig. 3 we plot variation 
of I with W for values of parameters as in Fig. 1. The 
plot exhibits a hell-shaped nature within a fairly narrow 
range of W around value 0.6. It may be noted that the 
minimum of a in Fig. 2 also occurs for W quite close to 
this value. 

We call the phenomenon observed above stochastic co- 
herence. This is similar to stochastic resonance phe- 
nomenon which shows a bell-shaped behaviour of tem- 
poral response as function of noise-strength Q . However 
one may note that our system does not have any intrinsic 
length-scale, whereas in stochastic resonance noise res- 
onates with a given time-scale; hence our use of word 
coherence rather than resonance [pi. In stochastic reso- 
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nance noise transfers energy to the system at a charac- 
teristic frequency, whereas in stochastic coherence noise 
induces coherence to the system. At this stage it is not 
clear to us as to how noise is inducing coherence to the 
system. However, absence of any length-scale over en- 
tire range of noise seems to indicate existence of weak 
self-organisation in the system induced by noise. 

In this context it should be noted that if a uniform- 
deviate random number is taken as site-variable, then the 
probability of a site to belong to a structure of length I is 
ps{l) « Z(2^)'~^(l — 2(5)^, where 5 is the structure param- 
eter introduced earlier. Hence number of such structures 
in a lattice of size L is given by n5(Z) « L(25)'~^(l — 2(5)^ . 
This shows an exponential decay of number of such struc- 
tures with length, which is to be contrasted with the 
power-law form (0) obtained for our lattice. 

It is possible to show how noise helps reducing instabil- 
ity of the structures, thereby increasing their abundance. 
Let us consider stability matrix M of a homogeneous 
state {■■■,Xt, Xt^Xt,---} (this may be thought of as a 
large structure with 5 = Q for simplicity). At time t -\- 1 
the matrix takes the simple form Mt+i — JFl , where 
J stands for the familiar tridiagonal matrix (with diago- 
nal elements 1 — e and offdiagonal elements e/2 on either 
sides) and F/ = F'{xt)i= dF/dx) = /i(l - 2xt) ■ Af- 
ter two timesteps, we get the stability matrix for three 
time steps as 5*3 = Mt+iMt+2Mt+3 = J^F^F^^^F^^^ = 
J3J;V[(^ _ 2Ft){l - 2^iFtil - Ft)} + 6m?7?(1 - m) - 
2r,t{l + - 6Ft + 6F2)} - ifii^f - 2r;t+i(l - 2Ft) + 
4774774+1] . For delta-correlated noise with uniform dis- 
tribution and zero mean (which is the case here), after 
averaging the above expression over noise-distribution 
one gets nonzero contribution due to noise only from the 
term quadratic in rj : 

< S3 >- J'(F;)V[i - 2fiFt{l - Ft) + e^i <772 > ] , (3) 

where <> denotes averaging over noise. For our lattice 
with localised structures in backdrop of spatiotemporal 
chaos we found the invariant density to be no longer sym- 
metric, with larger weightage for values of variable above 
0.5. Averaging expression ^ over this density makes the 
term 1 — 2/iFt(l — .Fi) negative. Adding positive contribu- 
tion of noise to it results in reduction of eigenvalue of the 
matrix S3 and hence consequent reduction in instability 
of the state. 

To study evolutionary aspects of these coherent struc- 
tures we obtained distribution of number n^r) of struc- 
tures vs. their lifetimes r for different W. This is shown 
in Fig. 4 (on log-linear scale). It exhibits decrease of n{T) 
with T with a stretched exponential type of decay having 
a form 



n(r) oc exp ( — (const. )r'') , 



(4) 



where /3 depends on W. We define average lifetime f of 
a structure as f = '''"-(''')/ S '^(''')- Fig. 5 we plot 
f vs. W for parameter values as in Fig. 4. The graph 



also shows a bell-shaped feature with maximum for W 
around 0.6. This W value is close to that corresponding 
to the extrema in figures 2 and 3. 

In order to ascertain chaotic nature of system evolu- 
tion we have calculated lyapunov exponent spectra for 
our system. We find a number of lyapunov exponents 
to be positive, implying that the underlying evolution 
is chaotic. Maximum lyapunov exponent Amax shows a 
minimum around noise-strength 0.6 ||]. We have stud- 
ied variation of A-spectrum with coupling strength e. 
Amax remains fairly constant for 0.2 < e < 0.8 for en- 
tire range of W. Behaviour of I with e is also investi- 
gated. It is found that I increases monotonically with 
e for all W. This fact is quite contradictory to what is 
expected from Amax- This indicates that lyapunov expo- 
nent alone cannot be used for proper characterisation of 
spatio-temporal features of the system. We have also ob- 
tained power spectrum of the time-series of the dynami- 
cal variable for a given site. The plot does not show any 
characteristic frequency and supports chaotic behaviour. 
In addition, we have calculated power spectrum of val- 
ues xt{i), i ~ 1, • • • , L, at a given time. Nature of the 
spectrum confirms our earlier observation that the sys- 
tem does not have any intrinsic length-scale. It may 
appear that the behaviour of our system is similar to 
spatio-temporal intermittency phenomenon llo[| . How- 
ever, our system does not show any regular burst-type 
feature (indicative of intermittency) in time; as already 
noted, power spectrum of time-series does not have any 
distinctive peak for entire range of noise-strength. Thus 
what we observe is development of appreciable coherence 
induced by noise in the system undergoing spatially in- 
termittent and temporally chaotic evolution. 

On quite a few occasions the entire lattice is seen to 
evolve as a single coherent structure to within the struc- 
ture parameter S. We have studied lifetime of these 
coherent structures. Full lattice coherence appears for 
W above a value around 0.4. Interesting thing is that 
occasionally this state is seen to persist for fairly long 
durations (at times as long as 200 timesteps or more), 
but eventually the coherence is seen to break up. This 
demonstrates that for our system synchronised state is 
not a stable attractor. 

However, in rare instances our system has been found 
to get into an apparently synchronised state after a very 
large time (larger than 10^ steps). This is obtained be- 
cause of finite accuracy of computation which cannot dis- 
tinguish unstable synchronised state from stable one . 
As demonstrated in the following, for our system this 
type of behaviour is an artifact resulting due to com- 
bination of finite size of the lattice and finite accuracy. 
Let us denote by T the average of timesteps required 
for first occurrence of full lattice coherence, the average 
taken over different initial conditions. We have stud- 
ied variation of T with structure parameter 5 for a fixed 
lattice-size L, which shows a definitive power-law of the 
form Ti (X S^'^ , where 7 is power-law exponent. This 
indicates that the synchronised state, i.e., full lattice co- 
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herent structure with ^ = 0, cannot be reached. We have 
also investigated T vs. L variation for a fixed 5. We see 
a power-law behaviour of form 

fsocL"" , (5) 
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where p is power-law exponent. Thus even full lattice 
coherence with finite 6 cannot be achieved as L ^ cx). 

We now show that existence of full-lattice coherence 
with nonzero S (distinct from synchronised state) in a 
finite lattice is essentially a consequence of power-law 
variation (^) of n{l) with /. From relation (^) proba- 
bility of a site to belong to a structure of length / is 
ps{l) oc = Therefore, the probability of a site 

to belong to a structure of length > L is Pq = Ps{> L) (x 

Y.T ~ ir °^ (fo'' " > 2, which is the 

case in the entire range of our observation as can be seen 
from Fig. 2). The probability that for the first time a 
site belongs to a structure of length > L at timestep T 
is then ^^(T) oc (1 - Po)^-iPo. Therefore average of T 
is given as fsi> L) = t.T=iT PsiT) oc (Po)"\ i.e., 

fs{> L) oc L"-2 . (6) 

This demonstrates that full-lattice coherence (with fi- 
nite 5) will be observed in a finite lattice. Fig. 2 tells 
us that a — 2 « 0.22 for W = 0.6, whereas we obtain 
i/(relation(|)) w 0.3. We beheve the discrepancy to be 
due to boundary corrections, since (^ is obtained for an 
infinite lattice whereas (||) holds for finite lattices. 

Above observations were also carried out for several 
values of e ranging from 0.1 to 0.9, as well as for nonlin- 
earity parameter fj, between 3.6 and 4. All the features 
remain essentially the same. We have also observed sim- 
ilar behaviour with periodic-boundary conditions for the 
lattice. 

In conclusion, we have reported a new phenomenon 
observed in a chaotically evolving one-dimensional CML 
driven by identical noise, which we termed stochastic co- 
herence. It is observed that there is a phenomenal in- 
crease in abundance of coherent structures of all scales 
due to noise. By considering stability matrix during three 
time steps, we have been able to show that noise can re- 
duce instability of these structures. Distribution of these 
structures shows a power-law decay with length of the 
structure, with an exponent which shows a minimum at 
some intermediate noise-strength. Average length as well 
as average lifetime of these structures exhibit character- 
istic maxima at a noise-strength quite close to previ- 
ous value. This bell-shaped feature is similar to that of 
stochastic resonance which is a temporal phenomenon. 
However, we emphasise that our system does not have 
any intrinsic length-scale, whereas stochastic resonance 
is associated with a particular time-scale. These obser- 
vations demonstrate that noise can play a major role in 
formation as well as in evolutionary dynamics of struc- 
tures in spatially extended systems. 
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FIGURE CAPTIONS 



Fig. 2. Variation of exponent a (relation (g)) with noise- 
strength W plotted for parameter values as in 
Fig. 1. 



Fig. 1. Plot of variation of number n{l) of structures with 
length I for a lattice with size L — 10000, for dif- 
ferent values of noise-strength W as indicated. Pa- ^ig 
rameters chosen are coupling strength parameter 
e = 0.6, structure parameter S = 0.0001, and non- 
linearity parameter /i = 4. Open-boundary condi- Fig. 
tions are used. Data are obtained for 100000 iter- 
ates per initial condition and 4 initial conditions. 



3. Plot of variation of average length I of structure 
with noise-strength W, with parameters as stated 
in Fig. 1. 

4. Plot of variation of number n(r) of structures with 
lifetime r, with conditions as in Fig. 1. 

5. Variation of average lifetime f of structures with 
noise-strength W shown for parameters as stated 
earlier. 
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Fig. 3 
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Fig. 5 
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